Nonlocally looking equations can make nonlinear quantum dynamics local 
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A general method for extending a non-dissipative nonlinear Schrodinger and Liouville-von Neu- 
mann 1-particle dynamics to an arbitrary number of particles is described. It is shown at a general 
level that the dynamics so obtained is completely separable, which is the strongest condition one can 
impose on dynamics of composite systems. It requires that for all initial states (entangled or not) a 
OO ' subsystem not only cannot be influenced by any action undertaken by an observer in a separated sys- 

0^ ' tem (strong separability), but additionally that the self-consistency condition Tr2 o0i+2 = 01 oTr2 

is fulfilled. It is shown that a correct extension to N particles involves integro-differential equations 
which, in spite of their nonlocal appearance, make the theory fully local. As a consequence a much 
larger class of nonlinearities satisfying the complete separability condition is allowed than has been 
assumed so far. In particular all nonlinearities of the form F{\il){x)\) are acceptable. This shows 
that the locality condition does not single out logarithmic or 1-homeogeneous nonlinearities. 
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I. INTRODUCTION 



> 

(N 

I In spite of the linearity of the Schrodinger and Liouville-von Neumann equations, nonlinear ly evolving states are 
. encountered in quantum mechanics quite often. Typically this is a result of approximations used in a description of 
collective phenomena (cf. nonlinear interferences in a Bose-Einstein condensate ) but no proof has been given so far 
that it is not the linearity that is a result of some approximation. This obvious observation motivated many authors to 
either look for nonlinear extensions of quantum mechanics, or to try to find an argument against a nonlinear evolution 
at a fundamental level. 

A privileged role in both cases was played by a separability condition. Apparently the first use of the condition can 
be found in The authors required that a nonlinearity must allow two separated, noninteracting and uncorrelated 
subsystems to evolve independently of each other, and this (plus some additional assumptions) led them to the 
^ nonlinear term In (jy'(x)|)?A(a;). Similar argument (with different additional assumptions) led Haag and Bannier 
and Weinberg [Q to the 1-homogeneity condition (i.e. ip{x) and X-tpi^x) should satisfy the same equation). An 
important element of these works was the assumption that the systems are uncorrelated, that is their wave function 
. ^ , is a product one. Equations that have such a locality property may be called weakly separable Q. It was shown 
^ ' later by Gisin and others that a weakly separable equation may still violate causality if an initial state of the 

5h composite system is entangled. A slightly different in spirit analysis of separability conditions was given in [ pl| where 
. . ' it was argued that an extension from 1 to particles may lead to new effects that become visible for N > Nq, where 
A^o is a parameter characterizing a given hierarchy of theories. 

It was originally conjectured that these properties refute any deterministic (i.e. non-stochastic) nonlinear gen- 
eralization of quantum mechanics. That this is not the case was shown for Weinberg-type quantum mechanics by 
Polchinski (for pure states) and Jordan (for density matrices). These results were subsequently generalized 
to a more general class of theories (Lie-Nambu dynamics) by myself in [p^ . Therefore, contrary to a rather com- 
mon belief, there exists a class of nonlinear generalizations of quantum mechanics that does not lead to the locality 
problems for general (pure-entangled and general-mixed) states. Such theories can be called strongly separable. The 
strongly separable theories considered so far involve equations which are Hamiltonian, which means there exists a 
Hamiltonian function that generates the dynamics. Actually, the solution of the problem proposed by Polchinski in 
[ p^ was based on an appropriate choice of this function. In addition, the Polchinski approach was formulated within 
a finite-dimensional framework. Still, there exists an interesting class of non- Hamiltonian and infinite-dimensional 
nonlinear equations. For example, it is known that all Doebner-Goldin equations that are Hamiltonian are lineariz- 
able ||l^ and therefore their strong separability may not be very interesting. ^Ouite recently the problem of strong 
separability of non-Hamiltonian Doebner-Goldin equations was addressed in [ |l6| . Using the typical extension of the 
dynamics from 1 to A^ particles the authors showed that at t = the time derivatives (up to the 3rd order) of a 
1-particle probability density in position space do not depend on the potential applied to the other particle if the 
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equation is Galilean-covariant. This suggests that Gahlean covariance may lead to strong locahty in nonrelativistic 
domain. The result agrees with computer algebraic tests undertaken by Werner who chose oscillator-type po- 
tentials and wave functions satisfying certain Gaussian Ansatz. The approach chosen in those works is essentially a 
third order perturbation theory applied to diagonal elements of reduced density matrices. A technical detail that did 
not allow these authors to find a general non-perturbative solution was that they investigated a dynamics of diagonal 
elements of a non-pure reduced density matrix (typical of entangled states) but the formalism they used was devised 
for state vectors and, hence, they had no control over the behavior of non-pure states and off-diagonal elements of 
density matrices |p^ . 

In this paper I will give a general and model independent solution of the problem. I will generalize to non- 
Hamiltonian systems the technique that proved useful in the context of density matrix equations, namely the triple- 



bracket formalism 21|. I will then show how to extend a 1-particle dynamics to A^-particle systems in a strongly 
separable way for a large class of both Hamiltonian and non-Hamiltonian Schrodinger non-dissipative equations. As 
opposed to basically all previous papers dealing with pure state nonlinear dynamics (the exception is , but it deals 
with finite dimensional systems) I will not guess the "obvious" form of the extension but will derive it. This somewhat 
long path beginning with pure states, then extension via mixtures, and again back to pure states, will be rewarded 
because the A^-particle form we will get will not be the one one might expect. Its paricularly interesting feature is 
the fact that the equations are integro- differential. This feature is implicitly present also in the Polchinski-Jordan 
formalism but is hidden behind the finite dimensional convention where the integrals do not explicitly show up. The 
approach is applicable to all non-dissipative Doebner-Goldin equations the equations discussed in H ^,22|, as 

well as a large class of nonlinear Schrodinger equations that were not discussed so far. 

I will prove simultanously a much stronger result. By the very construction the formalism is applicable to those 
Liouville-von Neumann nonlinear equations that reduce to the corresponding nonlinear Schrodinger dynamics on pure 
states. All these equations will be shown to be not only strongly separable but also completely separable. By the latter 
I mean a strongly separable dynamics which additionally satisfies the self-consistency condition 

Tr2 o0*_^2 = 01 °Tr2, (1) 

where 0*^2 ^^'^ 0i denote the dynamics of the composite system and the subsystem, respectively, and the partial 
trace Tr 2 is a map that reduces the dynamics from the large system to the subsystem. Condition (|l|) , which is 
independent of strong separability, was recently pointed out as an important ingredient of nonlinear dynamics that 
may be regarded as completely positive . We will also see that not only the probability densities in position space 
but also the off-diagonal elements of reduced density matrices are independedent of details of interaction in the remote 
systems. 



II. EXAMPLE: HAAG-BANNIER EQUATION AND ITS ALMOST-LIE-POISSON FORM 

The general convention I will use was elaborated in detail in but to make this work self-contained let us 

first explain the general scheme on a concrete example, the Haag-Bannier equation Q. This equation is of the 
Doebner-Goldin type, that is contains a 1-homogeneous nonlinear term with derivatives, but is simpler: 

^n^Ma) = (-I^A. + Via))Ha) + i(a) ^^^^1^^^^ (2) 
Its Liouville-von Neumann counterpart is 

ihdtp(a,a ) = -— Aa + V{a)]p[a,a ) + A[a)- -— r p[a,a) 

\ Zm / zipia, a) 

-[- — ^a' + V{a)jp{a,a) ~ A{a )^ p{a,a), (3) 

where p{a, a') — p{a', a), and /(a) = p(a, a) > is d"^a-integrable |2^ together with all its natural powers f{a)". For 
p{a, a') — xl){a)'ip{a') (||) reduces to the Schrodinger-Haag-Bannier dynamics. Using the composite index convention 
described in |14| we can write the equation in a form which is compact and simplifies general calculations. Denote 
Pa = p{a, a') and 

uat \ u! I \ iri ' \ ^ \n \x/ '\ 1 ^ /c^^y^(a-^y)Va[p(a,^/) ~ piy,a)] 

H (p) = H{a ,a) = Kia , a) + V(a)o{a — a) + A[a) — 7 oia — a). (4) 

zipia, a) 
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The kinetic kernel satisfies 

/ dyK{a,y)p{y,a') = - — Aap{a,a'), (5) 
J 2m 



and K{a,b) — K{b,a). All the integrals are in il^ i.e. da = d^a, etc., and the "summmation convention" is applied 
at the composite index level (two repeated indices are inegrated). There is no conflict of notation here because the 
composite indices are always in their upper or lower positions whereas the 3-dimcnsional coordinates a are arguments 
of functions or distributions. (Notice that the composite indices correspond always to pairs of primed and unprimed 
3-dimensional coordinates.) The indices are raised by the metric g"'' working as follows 

g^'^pf, = [ dbdb'S{a - b')S{b - a)p{b, b') = p{a, a) = 



So if Pa = p{a, a') then p" = p{a' , a). To lower an index one uses obvious inverse formul as p4| . 
The Liouville-von Neumann equation can be written in a triple-bracket-type form ]r^ , p6| , p7[ 

ihdtpf^ I dbdb'dcdc' (^5{f - b')6{b - c')S{c - f) - S{f - c')6{b - f')S{c - b')^ H{b' ,b) p{c' ,c) (6) 



HI' P" 



= %,ci?V = f^abc^ifV. (7) 
OPa 

Let us note that if 

H^^'Ji (8) 

opb 

then (0) describes a Lie-Poisson dynamics of a density matrix . If no Hamiltonian function H satisfying (H) exists, 
the dynamics will be called almost-Lie-Poisson. The generic case discussed in this Letter is almost-Lie-Poisson. The 
Weinberg-type and mean-field dynamics discussed in [p^p8| are Lie-Poisson. The discussion presented here applies 
to a general almost-Lie-Poisson dynamics where the structure constants have the following general form p^ ] 

^abc — -^af3^ -^f3j' -^ja' ^ a^' I fHa' -^^(5' (9) 
^abc ^ _^"/3'^/37'^7a' + ,^«7'^/3"'^7/3' (^Q) 

where uj°"^ = w° and laa' — la are, respectively, the symplectic form and the Poisson tensor corresponding to the 
pure state equation. Here uj"" — S{a — a'), la ~ S{a — a') but (^ and ( [lO| ) are valid also for other Hilbert spaces and 
equations (cf. 011). 



III. 2-PARTICLE EXTENSION 



Consider now a 2-particle system described by the density matrix which, depending on whether the state is pure 
or general, is in either of the two forms 

Pa = Paia2 = p{ai,a2,a[,a'2) (11) 

P=>(ai,a2)*(a'i,4). (12) 

The reduced density matrices are 

Pai^ j dy piai,y,a[,y) ^ u;"-''pa,a2 (13) 

pZ^ / d.ypiy,a2,y,a'2)=uj'''paia2- (14) 



The 2-particlc Liouville-von Neumann-Haag-Bannier equation is given by (0) but with all indices doubled like in ([Tl]), 
( p^ ) with the 2-particle structure constants given explicitly by 
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^i'l = Siai - b[)5{b, - c[)5{c, - a[) 
xS{a2 - b'2)5{b2 - C2)S{c2 - a'2) 
- S{ai - c[)S{bi - a[)S{ci ~ b[) 

x6{a2^c'2)S{b2-a'2)Sic2-b'2). (15) 

The crucial element of the whole construction is the nonlinear extension of the 1-particle Hamiltonian operator kernel. 
We define it as 

H\p)=Hf^ip^)p^ + Hfn/')^-^ (16) 
dpb dpb 

= H',^{p')u;'^+u;'-H',}{p''), (17) 

with b = &162. ( [T6| ) is a nonlinear functional generalization of the well known linear recipe -ffi+2 = Hi 1 + 1 ® and 
reduces to the formulas arising from the triple-bracket formalism if Hamiltonian functions exist. The whole 2-particle 
equation can be written as follows 

where P/ = /3/i/2 is the 2-particle density matrix. 

Now comes the important general result. Let us perform a partial trace i.e. contract both sides of ( p^ ) with lo^^ . 
Since pj_^ = i^^^ Pf^f^ — If^PfJ^ and — S{f2 — /2) is t- and p-independent 



ifidtp]^ 



OPa Opb 



KZi7^Hf^ip')i^p'' (20) 



fXcHf^ip')p'''- (21) 



where the transition from (|T^ to ( pO| ) is a consequence of 



aiN) ^P^d 5p" 

abc r- r- 

OPa Opb 



(22) 



holding for all A^-particle structure constants and reduced density matrices of non-overlapping subsystems (for the 
proof of (H^) see Lemma 1 in JT^). 

Let me summarize what has happened until now: Using a correct extension of a 1-particle nonlinear Hamiltonian 
operator and the general property of triple brackets we have reduced a 2-particle equation for a 2-particle density 
matrix to a 1-particle equation which involves only the quantities which are intrinsic to this subsystem. All elements 
depending on the other subsystem have simply vanished. We have obtained this by performing only one operation — 
the partial trace over the "external" subsystem. 

Therefore the reduced density matrix in / does not depend on details of interaction in the separated system // and 
the dynamics is strongly separable. Actually, we have simultaneously obtained more. Indeed, we do not assume that 
we take the partial trace at "t = 0" . Therefore we can trace out the external system at any time and the reduced 
dynamics is indistinguishable from a dynamics defined entirely in terms of the subsystem and starting from the initial 
condition pi(0) — p^(0) — Tr 2p(0) where p{Q) is the initial condition for the large system. It proves that the dynamics 
so constructed is completely separable. 

The result we have obtained is completely general and works for all Schrodinger equations whose nonlinear Hamilto- 
nian operator kernels allow to write them in terms of 1-particle density matrices. Putting it differently, the construction 
works correctly if a given Schrodinger equation allows for an extension to an almost-Lie-Poisson Liouville-von Neu- 
mann equation. The example of the Haag-Bannier equation served only as a means of focusing our attention and 
making the discussion less abstract. 

So how does the 2-particle equation look explicitly? Beginning again with the general formula we can obtain 
immediately its model independent form just by performing abstract operations on the composite indices. We get 
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Returning for the sake of completeness to the Haag-Bannier case we can write it as 
ihdtp{ai,a2, a'^, aj) 



-—(A,, + - - A„, ] + yi(ai) + 1/2(02) - ^^i(a'i) - 1^2(4) 



+Ai{ai) 



-^2(a2j 



/ dyi 5{ai - yi)VaJ/ dzp{ai,z,yi, z) - / dz p{yi, z,ai, z)] 
2i J dz p{ai, z, ai, z) 

J dy2S{a2 - V2)^a2 [/ dz p{z, Q2, z, 7/2) - j dz p{z, ^2, 02)] 
2i J dz p{z,a2,z,a2) 
^ - yi)^a[[J dz p{a[,z,yi,z) ~ J dz p{yi, z,a[, z)] 

^ 2i J dz p{a[, z,a'i, z) 

, ,J dy25{a'2 - y2)^a'[J dz p{z,a'2,z,y2) - J dz p{z,y2, z,a'2)] 



2i J dz p{z,a2, z,a2) 



p(ai,a2,ai,a2). 



(23) 



Its pure state 2-particle counterpart is 



ihdt'^{ai,a2) 



2m 



(A„, +A,,) +Vi(ai) + y2(a2) 



+Ai{ai) 



+^2(02) 



/ dz[*(ai,z)Vai^'(ai,2) - *(ai, 2;)Vai^'(ai, z) 



2i /rfz|*(ai,z)|2 



/ dz [*(z, a2)Va2^'(z, 02) - ^(z, a2)Va2*(z, 02) 



2i / dz|^'(z,a2)p 



^-(01, 02). 



(24) 



This equation has several interesting features. The Hamiltonian operator obviously reduces to the sum of ordinary 
1-particle terms, involving no integrals, if 5* is a product state. What makes it unusual is the presence of the integrals. 
Typically it is said that such equations should not be taken into account because they are nonlocal. They indeed 
appear nonlocal but a closer look shows that it is in fact just the opposite: The currents and probability densities 
are the local 1-particle ones. Therefore it is the lack of appropriate inegrals that makes typical 2-particle equations 
nonlocal. The equation considered in ||l^ involved no such integrals and this led to difficulties. An extension of the 
above results from 2 to particles is immediate so explicit formulas will not be given. 



IV. FURTHER EXAMPLES 



Let me now list some of nonlinear Hamiltonian operators that have been considered in literature and which admit 
the completely separable extension to an arbitrary number of particles in arbitrary entangled and mixed states. The 
general rule is that an extension to N particles will be given in one of these forms but instead of p{x, x') a reduced 
1-particle density matrix should be placed. If the A^-particle state is pure then the reduced density matrix is a 
functional of the pure state which involves A^ — 1 integrations. This matrix is subsequently put into a suitable place in 
the Schrodinger equation. Similarly it can be put into a Liouville-von Neumann equation if the state is more general. 

a) "Nonlinear Schrodinger" 

|V'(a;)p p{x,x) 

b) Biahjnicki-Birula-Mycielski 

ln(|Vj(a;)p) -^\np{x,x) 
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Obviously in the same way one can treat any equation with nonhnearities given by some function F{\^p{x)\) [ pO[ 
c) Doebner-Goldin 



R 



1 • 



R2 



R^ 



R 



4 • 



R 



5 • 



1 ll){x)/S.x'4'{x) - 'lp{x)A3:'lp{x) 

1 J dz6{x ~ z)Ax[pix, z) - p{z,x)] 
2i p{x,x) ' 

Aa:|V'(a:)P _^ A^p{x,x) 
\ip{x)\'^ p(x,x) 

1 [Vi(x)V^^(x) - ^(x)V,^(x)] ' 

1 (/dz(5(a; - z)V2:[p(a;, z) - p(z,a;)])^ 
^ p(a;,a;)2 

1 / dz6{x - z)Va; [/9(a:, z) - p(z, a;)] • Vxp{x, x) 
2i p(a;,a;)2 

\ip{x)\'^ p{x,xY 



d) Twarock on 5^ p| 



ip{x)ip{xy — 4'{x)ip(x)' 

_^ [/ dy5{x - y)dlp{x, y)] [/ dz5{x - z)dlp{z, x)] - c.c. 
p{x, x) J dy6{x - y)dxp{x, y) - c.c. 

e) (n, n)-homogeneous nonlinearities . Denote by Z? a differential operator involving arbitrary mixed partial derivatives 
up to order k. Consider a real function F{ijj) = F(^Difi{x)) , (n, n)-homogeneous i.e. satisfying F{\ip) — X'^X"-F{ip). 
We first write 



F{^{x)D^{x)) 

and then apply the tricks used for the Haag-Bannier, Doebner-Goldin and Twarock terms. Obviously any reasonable 
function of such (n, n)-homogeneous expressions with different n's is acceptable as well. 

V. EXTENSION IN STAGES 

The formalism proposed above allows to extend a 1-particle dynamics to N particles. We will now show that it 
allows for a more general kind of extension: From 1 system to N systems. The procedure is self-consistent in the sense 
that one can first produce several composite systems from the 1-particle ones, and then combine them into a single 
overall composite system. Alternatively, one can produce the final system without the intermediate stages, directly 
by extension from single particles. 

Consider the A^-particle extension 

H'Ap)-H'j-''^ip) 

= H^iM'')r^ + Ht^M'')T^ + ■■■ + HlSM''^)P^ (25) 
= . . .cj''" + w^ii/['2yp(2)),^63 . . + . . . + cj''! . ..uj^«-^H\^^{p^^^). (26) 



6 



The Hamiltonian operator (kernel) i?^ (p) describes a composite system, labelled A, which consists of N particles that 
do not interact with one another. Consider now another system, labelled B, consisting of M particles. Its Hamiltonian 
operator is 

g {N+1) g {N+2) 

= Hf^l\y^+'^) '""-^ + Hl-XlM^^'^) + . . . + Hl^XtiM''^''^) '""^^ (27) 

oPbjv+i-.-biv + Af "Pb]V + l---bjv + M ''PbN+l-..6jv+M 

+ . . . + . . . cj''"+^'-ii?J';^+*/^)(p(^+*^)) (28) 

The {N + Af)-particle extension can be obtained in two stages 

X„A g„B 
TTdi...dM ( „A\ "(ii...(iiv I rrdjv + i---djv + Af / B\ '^<ijv + i---<i]v + Af 

"Hbi-.-bN + M ''Ffci ...bjv + M 



7V4-M 



= (H(\yp(i')cj''^ . . .cj^" + . . . + c^''^ . . w^"-ii/['^j(pW))cj''"+i . . . 

+uj^^ ...uj'"' (i/[';^+\^(p(^+i))w^"+^ . . . tj''"+" + . . . + . . . cj''"+'-^-ii7j';^+^^)(p(^+*-f))) (29) 

which could be derived also directly from the (N+M)-particle extension of the 1-particle dynamics. 



VI. IS NONLINEAR QUANTUM MECHANICS NONLOCAL? 

N. Gisin in his nowadays classic paper argued that any nonlinear and non-stochastic dynamics necessarily leads 
to faster-than-light communication via EPR-type correlations. The argument is based (implicitly) on the additional 
assumption that a reduced density matrix (of, say, Alice) evolves by means of an independent dynamics of its pure- 
state components. (Let us note that this was essentially the main question discussed by Haag and Bannicr in the 
context of the nonlinear, convex scheme proposed by Mielnik [^.) This viewpoint is suggested by the fact that what 
one starts with is a nonlinear (Schrodinger-type) dynamics of pure states. But in the EPR case the pure state is a 
2-particle one and one has to start with a 2-particle nonlinear Schrodinger equation. Assume this equation is given 
(e.g. our Haag-Bannier equation) and its 2-particle pure-state solution \'ip) has been found. To have the Gisin effect 
one has to find at least one observable of the form 

{A) = {'^\A® IIV-) = TtapaA (30) 

describing the Alice subsystem, which depends on a parameter controlled by Bob. Here pA = Trs|'0)(V'l is the 
reduced density matrix of the Alice subsystem. Assume that the only actions Bob can undertake are reducible to the 
modifications of the parameters of the Hamiltonian corresponding to "his" particle. (In the case contemplated by Gisin 
Bob would rotate his Stern-Gerlach device; this is equivalent to modifying the magnetic term in the corresponding 
interaction Hamiltonian.) (|3^) shows that Bob can influence {A) if and only if the reduced density matrix pA depends 
on the parameter he controls. This is true for the Weinberg theory [lO| ] but will never happen for any equation we 
have discussed. It follows that the Gisin effect applies to a limited class of theories. 

The additional assumption that leads to the Gisin phenomenon is the following: Each time Bob makes a measure- 
ment, he chooses an initial condition for the reduced density matrix of Alice. This explicitly involves the ordinary 
reduction of the wave packet postulate. One can invent other versions of the postulate, for example, a composition 
of a nonlinear map N with a projector P: o P o N, as proposed by Liicke in his analysis of nonlinear gauge 

transformations [Q. 

Alternatively, as in the Weinberg theory, one can get the Gisin effect without the explicit use of the projection 
postulate: The nonlinear generator of evolution must be basis dependent. If this is the case, the 2-particle solution 
is parametrized by the choice of the basis in the Hilbert space. Assuming that Bob can change this basis we allow 
him to change globally the form of the solution. But the rotation of the Stern-Gerlach device (modification of the 
infarction term) is not yet a change of basis in this sense. In the Weinberg theory the generator was interpreted as 
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an average energy. A change of basis was equivalent to a different way of measuring the average. For this reason it 
was justified to say ||l2| that the effect was imphed by the projection postulate although explicit calculations were 
not referring to any projections [ p^ (similarly to the linear case, the projection is not included in a non-stochastic 
Schrodinger dynamics). 

We can conclude that the Gisin-type reasoning has to be regarded as "unphysical" , exactly in the same sense as 
the Einstein-Podolsky-Rosen argument was unphysical according to Bohr. The version of nonlinear theory we have 
discussed seems to be free of physical inconsistencies and is not more nonlocal than the linear one. 

The idea of using the Haag-Bannier equation as a laboratory for testing the concepts of separbility was suggested 
to me by G. A. Goldin who also informed me about the results of his joint work with G. Svetlichny. I am indebted 
to G. A. Goldin, H.-D. Doebner, W. Liicke, P. Nattermann, and J. Hennig for valuable suggestions and critical 
comments, and to W. Puszkarz who explained to me how to include into the framework the equations of Kostin p3] 



and Staruszkiewicz |34 3q]. The work is a part of the Polish-Flemish project 007 and was done during my stay at the 



Arnold Sommerfeld Institute in Clausthal. I gratefully acknowledge a support from DA AD. 
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